Magnetic properties of spin diluted iron pnictides from /iSR and NMR in 

LaFei ^Ru^AsO 
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The effect of isoelectronic substitutions on the microscopic properties of LaFei-j^Ru^^AsO, for 
< a; < 0.8, has been investigated by means of fiSR and ^^^La NMR. It was found that Ru 
substitution causes a progressive reduction of the Neel temperature (Tjv) and of the magnetic order 
parameter without leading to the onset of superconductivity. The temperature dependence of ^'^^La 
nuclear spin- lattice relaxation rate l/Ti can be suitably described within a two-band model. One 
band giving rise to the spin density wave ground-state, while the other one is characterized by 
weakly correlated electrons. Fe for Ru substitution yields to a progressive decrease of the density of 
states at the Fermi level close to the one derived from band structure calculations. The reduction of 
Tn with doping follows the predictions of the Ji — J2 model on a square lattice, which appears to be 
an effective framework to describe the magnetic properties of the spin density wave ground-state. 

PACS numbers: 76.75.-|-i, 75.10.Jm, 74.90.-|-n, 76.60.Es 



I. INTRODUCTION 



The recent discovery of high temperature supercon- 
ductivity nearby the disruption of magnetic order in Fe- 
based compounds^ has stimulated the scientific commu- 
nity to further consider the role of magnetic excitations^ 
as a possible candidate for the pairing mechanism. In 
order to address this point an appropriate modellization 
of the microscopic magnetic properties of the Fe based 
superconductors and of their parent compounds is nec- 
essary. In the phase diagram of iron pnictides supercon- 
ductivity appears to compete with a commensurate spin 
density wave (SDW) magnetic order'' characterized by a 
reduced magnetic moment The magnitude of this mo- 
ment is much lower than that evaluated from band struc- 
ture calculations)^ possibly due either to the strong elec- 
tronic correlations, not appropriately taken into account 
in those calculations, and/or from frustration effects.^ 

When the SDW phase is suppressed, by chemical sub- 
stitution or by applying high pressures, superconductiv- 
ity is usually recovered i=2'"=§ In fact, in A(Fei_2;Ru2;)2As2 
(A=Sr, Ba) the substitution of Fc by Ru suppresses the 
magnetic ordering for x — > 0.3 and leads to bulk super- 
conductivity for 0.2 < a; < OA^^:^ On the other hand, it 
has been shown that in PrFei_a;Rua; AsOfi^ in spite of an 
analogous disruption of the SDW ordering for x ~ 0.6, no 
superconductivity is found up to the full Ru substitution. 
Indeed band structure calculations for REFei_a;Ru2;AsO 
show that the electronic structure around the Fermi level 
is actually only slightly affected by Ru substitution^- 
and that only a minor charge doping takes place, even 
in the presence of very large Ru contents. This means 
that in this system the Fe/Ru substitution is effectively 
isoelectronic and accordingly no relevant modification of 
the electronic ground-state is observed. Furthermore it 



has been predictedi^ that the Ru atoms do not sustain 
any magnetic moment suggesting that REFci-^jRuajAsO 
should be considered as a spin-diluted system. 

In order to understand which is the effect of Ru 
substitution in the 1111 family of iron pnictides we 
have performed yuSR and -'^^^La NMR measurements in 
LaFci_a;Rua;AsO. A progressive reduction of the Neel 
temperature (Tat) and of the magnetic order parameter 
is observed with increasing x. Both quantities eventu- 
ally vanish for x — 0.6 without leading to the onset of 
superconductivity, at least up to a; = 0.8. The tempera- 
ture dependence of ^^^La nuclear spin-lattice relaxation 
rate 1/Ti can be suitably described within a two-band 
model, one giving rise to the SDW ground-state, while 
the other one characterized by a Fermi-gas behaviour. 
Ru for Fe substitution yields to a progressive decrease of 
the density of states at the Fermi level, a trend which is 
quite consistent with band structure calculations. The 
low temperature behaviour of 1 /Ti in the ordered phase 
and the reduction of T/v with doping can both be de- 
scribed within the Ji — J2 model on a square lattice and 
indicate that LaFei_;rRu;rAsO behaves as a spin-diluted 
system with competing exchange interactions. This ob- 
servation suggests that this model provides an effective 
framework to suitably describe the role of frustration in 
the iron pnictides i^ii^ 



II. TECHNICAL ASPECTS AND 
EXPERIMENTAL RESULTS 

The experiments were performed on polycrystalline 
LaFei_a;Rua;AsO (0 < x < 0.8) samples prepared as de- 
scribed in Ref il7i. Structural characterization was per- 
formed by X-ray powder diffraction at room tempera- 
ture and Rietveld refinement was carried out on selected 
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diffraction patterns. Microstructure was inspected by 
scanning electron microscopy. Transport measurements 
sfiow a decrease in the resistivity p with, increasing Ru 
content and a shift to low temperature of the characteris- 
tic anomaly in dp/dT at T/v The magnetic susceptibil- 
ity X derived from SQUID magnetization measurements 
showed an analogous behaviour of the peak in dx/dT 
at T/v, shifting to lower temperatures with increasing x. 
However, for x > 0.4 the peak progressively smeared out 
and eventually for x — >■ 0.6 the susceptibility showed ba- 
sically a Curie- Weiss behaviour. 

Zero field (ZF) /iSR experiments have been performed 
at PSI with GPS spectrometer. For T > Tat, the muon 
asymmetry is characterized by a decay which progres- 
sively increases with decreasing temperature. Below Tjv, 
at low Ru contents, well defined oscillations are observed 
(FiglT]), evidencing the presence of a magnetic order. For 
X > 0.3 (FigH]) these oscillations are markedly damped 
due to the increase in the local field distribution at the 
muon site. 




FIG. 1: Time evolution of the muon asymmetry in 
LaFei_a;Rua;AsO for x = 0.1 (top) and for x = 0.4 (bottom) 
at three selected temperatures. The solid lines are the best 
fits according to Eq[l] 

Accordingly, for x < 0.3 the asymmetry could be fitted 
to the sum of a fast and a slow oscillation plus a non- 
oscillating term, namely 

A{t) = A,e-^^'fi-f^B^t)+A2e~^^'f{^^B^t)+A^^e-^ii' , 

(1) 

where 7^ is the muon gyromagnetic ratio, Bf^ 2 is the 
local field at the muon sitesi^ii^ « = 1 or 2 and A1.2 



are the corresponding decay rates. As can be seen from 
FigH] the amplitude of the fast oscillating component is 
significantly larger than the one of the slow oscillating 
term. The fast oscillating signal could be reproduced by 
fil^B'^t) = cos{-/f,B'^t + (f)) for x < 0.3. At higher Ru 
contents f{jf^B'^t) — Jo(7^i3ft), with Jq the zeroth or- 
der Bessel function. For T <^ Tn, for x < 0.6, one finds 
that A|| is 1/3 of the total asymmetry, as expected for 
fully magnetically ordered powdersi^ Then it is possible 
to estimate the temperature dependence of the magnetic 
volume fraction vm{T) — (3/2)(l — A||(T)) from the tem- 
perature dependence of A\\{T). The fraction vm{T) is 
shown in Fig. [2] for different doping levels. One notices 
a fast drop of vm for T — > T/v, which can be empirically 
fitted with vm{T) = 0.5(l-er/(T-r^7%/2Ay)), where 
T^" represents an average transition temperature. One 
observes a progressive decrease of with x together 
with some broadening, likely due to inhomogeneity in 
the Ru distribution. Notice that for a; < 0.5 the temper- 
ature at which vm reaches about 100 % is close to Tm as 
determined from resistivity measurements, as it is shown 
in details in the next section. 

The fit of the muon asymmetry for a; < 0.3 allows one 
to derive the temperature dependence of the modulus of 
the local field at the muon Bf (T) = A; < S >, with A 
the hyperfine coupling tensor and < S > the average Fe 
spin value, corresponding to the order parameter. In Fig. 
[3]the temperature dependence of i?f (T), the internal field 
at the muon site close to FeAs layersjii^ is reported for 
X < 0.3. One notices a progressive decrease of the Fe 
magnetic moment with increasing x. 
120- 
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FIG. 2: Temperature dependence of the magnetic volume 
fraction in LaFei_a;Rua; AsO. The solid lines are the best 
fits according to the phenomenological expression vm{T) = 
0.5(1 -er/(r-n7^A)). 

^•^^La NMR measurements have been carried out by us- 
ing standard radiofrequency pulse sequences. The spec- 
tra of the X — compound was in excellent agree- 
ment with the one previously reported by Ishida et al^ 
Upon doping one observes that the peaks of the cen- 
tral line progressively smear out and the spectrum gets 
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FIG. 3: Temperature dependence of the local field at the 
muon site f in LaFei-^^Rui^AsO, for 0.1 < x < 0.3. 
The dashed lines are guides to the eye. In the inset the 
a;— dependence of the T ~^ order parameter < 5 > (T — >■ 
0,x)/ < S > {T ^ 0,0) = Bf(r ^ 0,x)/B'^{T 0,0) is 
shown. The dashed line is a guide to the eye. 



narrower, suggesting a decrease in the electric field gra- 
dient at ^'^^La nuclei and possibly also a change in the 
paramagnetic shift tensor. The explanation of this phe- 
nomenology, however, goes beyond the aim of the present 
manuscript. Nuclear spin-lattice relaxation rate 1 /Ti was 
measured on the central transition by using a satura- 
tion recovery pulse sequence. The recovery of the nuclear 
magnetization does not follow the trend expected for a 
magnetic relaxation mechanismiS^ by assuming a single Ti 
but could rather be fitted by assuming a distribution of 
Ti. Hereafter Ti will be defined as the time at which the 
recovery law decays to 1/e. The temperature dependence 
of 1 /Ti in LaFei_a;Rua; AsO is shown in FiglH 1 /Ti shows 
a linear temperature dependence above Tat and a drop 
below. The peak at the transition temperature originates 
from the divergence of critical fluctuations, progressively 
smeared out by inhomogeneities for increasing x. Figure 
[S] displays the temperature dependence of l/TiT on log- 
linear scales. It is worth noticing that the relaxation rate 
does not vanish for T — 0, where it tends to a constant, 
cc— dependent value. The low temperature levelling of 
1/Ti has already been observed by Nakai et al. in Refl23l 
for the parent x = compound. 



III. ANALYSIS AND DISCUSSION 

First we shall concentrate on the effect of Ru doping on 
the sublattice magnetization and on the density of states 
at the Fermi level in the light of band structure calcu- 
lations, performed within the density functional theory 
(DFT) using the local density approximation (LDA) for 
the exchange and correlation functional. The band struc- 
ture was obtained with the Siesta code,^** which utilizes 
a linear combination of atomic orbitals for valence elec- 
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FIG. 4: Temperature dependence of ^^^La 1/Ti in 
LaFei_a;Rua;AsO, for 0.1 < a; < 0.6, in a 7 Tesla magnetic 
field. 
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FIG. 5: Temperature dependence of "^La 1/TiT in 
LaFei_iRu2:AsO, for 0.1 < a; < 0.6, in a 7 Tesla magnetic 
field. 



trons and separable norm conserving pseudopotentials 
with partial corrections for atomic cores. Crystal struc- 
tures relaxations were performed with periodic boundary 
conditions. The atomic positions as well as the cell struc- 
ture were allowed to be optimized in the paramagnetic 
state by using a conjugate gradient algorithm. The real 
space integration grid had a cut-off of 500 Ry and up to 
12000 points were used for the Brillouin zone sampling 
using the Monkhorst-Pack k-points sampling. Stringent 
criteria were adopted for the electronic structure con- 
vergence and equilibrium geometry (residual forces lower 
than 10-2 eV/A). 

Our results are in close agreement with previous band 
structure calculations reported in Ref. [HI. We find that 
Fe magnetic moment decreases with Ru doping and even- 
tually vanishes for x around 0.5. This trend is quali- 
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tatively consistent with the one experimentaUy derived 
from the x dependence of Bf(r 0,x) (FigE]). Nev- 
ertheless, as it has already been pointed out in the in- 
troduction, band structure calculations do not provide 
a quantitatively correct estimate of the magnitude of Fe 
magnetic moment, as well as of the a;-dependence of the 
order parameter. Namely, the initial slope d < S{T — 
0, x) > / dx obtained experimentally is much faster than 
the one derived from band structure calculations. 

On the other hand, a better agreement with the ex- 
perimental findings is observed for the calculated density 
of states at the Fermi level. This quantity can be de- 
rived experimentally from 1/Ti measurements 3^ In fact, 
above Tn l/Ti follows the Korringa behaviour expected 
for a Fermi liquid (Fig]?]), namely 1/Ti = CuqT, with hq 
the density of states at the Fermi level and C a constant 
accounting for the hyperfine coupling between the elec- 
trons and ^'^^La nuclei. Then, by taking the value of the 
spin-lattice relaxation rate around room temperature one 
can write that y^Ti{0)/Ti{x) ~ no(a;)/no(0) and derive 
the x-dependence of the density of states at the Fermi 
level. In FiglB] this ratio is compared to the results ob- 
tained from ab initio calculations in the paramagnetic 
state. One observes that no (a;) decreases with increasing 
Ru content, in excellent agreement with band structure 
calculations. This decrease should be associated with the 
larger extension of Ru d orbitals which leads to an en- 
hanced delocalization of the electrons. 

At low temperature, in the SDW phase, LaFeAsO nu- 
clear spin-lattice relaxation can be described as the sum 
of two contributions, as suggested by Smerald et al^ 
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where l/Tf ^ cx T is a Fermi-gas like term arising from 
weakly correlated electrons which accounts for the low- 
temperature levelling of l/TiT, while l/Tf'^ is the one 
from a band with strongly correlated electrons where the 
Fermi surface nesting leads to the insurgence of the SDW 
phase. 

At low temperature l/Tf^ reduction is determined 
by the gap A in the spin excitations. In fact, following 
Refl26l. taking into account that the experiments were 
performed on powders, one can write 
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where ^[x] 



KbT 



(3) 



:2Lii(e-i/")+a;3Lz2(e-i/") 



KbT 



with Lin{z) the n polylogarithm of z. In Eql3]y^h is 
the hyperfine coupling between the longitudinal fluctu- 
ations of the Fe moment and the nuclear spin, 77V is 
"'^'^^La gyromagnetic ratio, toq the amplitude of the fluc- 
tuating moment, while a = V/x±v^''^^^'^ , with V the 
unit cell volume, Vs the average spin-wave velocity and 




FIG. 6: The density of states at the Fermi level is reported 
as a function of Ru content in LaFei-ajRu^^AsO, normalized 
to the a: = value. The squares represent the data derived 
from '^^^La Ti measurements, while the circles are from band 
structure calculations. The dashed line is a guide to the eye. 



X± the transverse spin susceptibility. It is interesting to 
observe that this approach applies also to the Ji — J2 
model on a square lattice with localized spins.— In fact, 
although LaFei_a;Rua;AsO is not a localized spin system, 
the Ji — J2 model appears to be still applicable in some 
effective version also to the iron pnictides. Further sup- 
port to this idea will be presented subsequently in the 
discussion of LaFei^ajRu^jAsO phase diagram. 

It is interesting to notice that many parameters ap- 
pearing in Eq|3] determine also the reduction of the sub- 
lattice magnetization. In fact one has that 



< 5" > (0)- <S>{T) 
< 5 > (0) 



3/2 



-A/ksT 



(4) 

Thus, both the temperature dependence of < S" > (T) 
and of 1/Ti for T ^ are determined by A. By fitting 
ZF-^SR and ^^^La 1/Ti curves for x = 0.1 (FigEl) one 
finds an accurate fit of both quantities for A = 8 ± 2 
meV, a value which is close to the one found in other iron 
pnictides by means of inelastic neutron scattering 1^"— 
By assuming in Eql3] mo ~ l/is one finds a quantitative 
agreement with ^'^^La 1/Ti data for Ah — 1-8 kG/^B- 
This value is reasonably close to that estimated in the 
SDW phase from the hyperfine field at the ^^^La nucleus 
B^a — 2.4 kG. Namely, by assuming the T = value 
of the magnetic moment of Fe in the SDW phase'^^ to 
be m = 0.6^5, one has Ah — \BLa\lz < M(0) >~ 1 
kG//iB, with z = 4 the nearest neighbor Fe atoms. 

Further insights on the applicability of the J1-J2 model 
to LaFei_a;Rua;AsO comes from the analysis of the phase 
diagram. As it is shown in Fig[8] the Neel temperature 
determined either from transport or from ZF-/iSR by tak- 
ing the temperature at which vm — >■ 1 (Fig. 2), decreases 
almost linearly with x and eventually vanishes around 
X = 0.6. It may be argued that x ^ samples present 
a distribution of Neel temperatures and that the criteria 
chosen for identifying Tn may vary with the determina- 
tion technique. Still, if we choose, for example, to evalu- 
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FIG. 7: (Top) Temperature dependence of ^^^La 1/Ti in 
LaFeo.gRuo.iAsO, in a 7 Tesla magnetic field, for T <^ Tn ■ 
The solid line shows the best fit according to Eq|3] (Bottom) 
Temperature dependence of the magnetic order parameter de- 
rived from ZF-/iSR measurements. The solid line is the best 
fit according to Eq|4l 



model on a square lattice. 



IV. CONCLUSIONS 

In conclusion, we have shown that Ru for Fe substitu- 
tion in LaFei_a;Rua;AsO causes a progressive reduction of 
the Neel temperature (Tjv) and of the magnetic order pa- 




FIG. 8: The doping dependence of Neel temperature in 
LaFei-xRuiiAsO, normalized to its a; = value is shown. Cir- 
cles show the behaviour derived from transport data, squares 
from ZF-/jSR, while diamonds show the corresponding be- 
haviour for Li2Vi_a;Tia;Si05. 



ate the average T^" from the flex in the muon magnetic 
fraction Vm{T), the slope dT^^/dx turns out to be the 
same as that of Fig. [51 This points out that our main 
findings are not affected by some distribution in the Neel 
temperatures. 

Now, it is instructive to compare the slope dTM{x)/dx 
with the one found in Li2Vi_a;Tij:Si05, a prototype of the 
Ji — J2 model on a square lattice with spin dilution arising 
from the substitution of V'^+ {S = 1/2) with {S = 

0). One notices that the trend of TN{x)/Tpf{0) is the 
same in LaFei-^^Ru^jAsO and Li2Vi_2:Tia;Si05, further 
supporting the idea that the Ji — J2 model on a square 
lattice^^ provides an effective framework to appropriately 
describe the effects of the spin dilution induced by Ru 
for Fe substitution.^^ Moreover, it is pointed out that in 
LaFei_2:Ruj;AsO Tm{x) vanishes for x around 0.6 which 
corresponds to the percolation threshold for the Ji — J2 



rameter without leading to the onset of superconductiv- 
ity. The analysis of ^^^^La nuclear spin-lattice relaxation 
rate 1/Ti indicates that this system can be described 
within a two-band model, one of them giving rise to the 
spin density wave (SDW) ground-state. Fe for Ru sub- 
stitution yields to a progressive decrease of the density of 
states at the Fermi level in quantitative agreement with 
band structure calculations. The behaviour of 1/Ti in 
the SDW phase and the reduction of T/v with Ru substi- 
tution can both be described within the Ji — J2 model on 
a square lattice which suggests that LaFei_a:Rua:AsO be- 
haves as a spin-diluted system with competing exchange 
interactions, pointing out the relevant role of frustration 
in the parent compounds of iron based superconductors. 
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surements at PSI and the access support by EU con- 
tract RII3-CT-2003-505925 (NMI3) are gratefully ac- 
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